function [VecResults,VecOutliers2]=RemoveOutliers(VecTemp,ValueTemp4)

VecOutliers2=[];
ValueMean=mean(VecTemp);
ValueSQRT=sqrt(var(VecTemp));
ValueTemp=find(VecTemp<ValueMean+ValueTemp4*ValueSQRT&VecTemp>ValueMean-ValueTemp4*ValueSQRT);
VecOutliersIndexTemp=ValueTemp;
ValueTemp3=find(VecTemp>=ValueMean+ValueTemp4*ValueSQRT);
VecOutliers2=[VecOutliers2;VecTemp(ValueTemp3,:)];
ValueTemp3=find(VecTemp<=ValueMean-ValueTemp4*ValueSQRT|VecTemp>=ValueMean+ValueTemp4*ValueSQRT);
VecOutliersIndex=ValueTemp3;
VecTemp=VecTemp(ValueTemp,:);
while size(ValueTemp3,1)>=1
    ValueMean=nanmean(VecTemp);
    ValueSQRT=sqrt(nanvar(VecTemp));
    ValueTemp=find(VecTemp<ValueMean+ValueTemp4*ValueSQRT&VecTemp>ValueMean-ValueTemp4*ValueSQRT);
    ValueTemp3=find(VecTemp>=ValueMean+ValueTemp4*ValueSQRT);
    VecOutliers2=[VecOutliers2;VecTemp(ValueTemp3,:)];
    ValueTemp3=find(VecTemp<=ValueMean-ValueTemp4*ValueSQRT|VecTemp>=ValueMean+ValueTemp4*ValueSQRT);
%     VecOutliersIndex=[VecOutliersIndex;VecOutliersIndexTemp(ValueTemp3,1)];
%     VecOutliersIndexTemp=VecOutliersIndexTemp(ValueTemp,:);
    VecTemp=VecTemp(ValueTemp,:);
end
VecResults=VecTemp;